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Abstract 

Poulalhon and Schaeffer introduced an elegant method to linearly encode a pla¬ 
nar triangulation optimally. The method is based on performing a special depth-first 
search algorithm on a particular orientation of the triangulation: the minimal Schny- 
der wood. Recent progress toward generalizing Schnyder woods to higher genus 
enables us to generalize this method to the toroidal case. In the plane, the method 
leads to a bijection between planar triangulations and some particular trees. For 
the torus we obtain a similar bijection but with particular unicellular maps (maps 
with only one face). 


1 Introduction 

A closed curve on a surface is contractible if it can be continuously transformed into a 
single point. In this paper, we consider graphs embedded on a surface such that they 
do not have contractible cycles of size 1 or 2 (i.e. no contractible loops and no multiple 
edges forming a contractible cycle). Note that this is a weaker assumption, than the 
graph being simple, i.e. not having any cycle of size 1 or 2 (i.e. no loops and no multiple 
edges). A graph embedded on a surface is called a map on this surface if all its faces are 
homeomorphic to open disks. A map is a triangulation if all its faces are triangles. We 
denote by n the number of vertices, m the number of edges and / the number of faces 
of a given map. 

Poulalhon and Schaeffer introduced in [18] a method (called here PS method for 
short) to linearly encode a planar triangulation with a binary word of length log 2 (^) ~ 

*This work was supported by the grant EGOS ANR-12-JS02-002-01 and the project-team GALOIS 
supported by LabEx PERSYVAL-Lab ANR-ll-LABX-0025. 

^ Vincent .despre@gipsa-lab .fr 

Ganiel.goncalves@lirmm.fr 

Genjamin.leveque@lirmm.fr 


1 





n log 2 (^) ~ 3, 2451 n bits. This is asymptotically optimal since it matches the infor¬ 
mation theory lower bound. The method is the following. Given a planar triangulation 
G, it considers the minimal Schnyder wood of G (that is the orientation where all inner 
vertices have outdegree 3 and that contains no cycle oriented clockwise). Then a special 
depth-first search algorithm is applied by “following” ingoing edges and “cutting” out¬ 
going ones. The algorithm outputs a rooted spanning tree with exactly two leaves (also 
called stems) on each vertex from which the original triangulation can be recovered in 
a straightforward way. This tree can be encoded very efficiently. A nice aspect of this 
work, besides its interesting encoding properties, is that the method gives a bijection 
between planar triangulations and a particular type of plane trees. 

Castelli Aleardi, Fusy and Lewiner [5] adapt PS method to encode planar triangula¬ 
tions with boundaries. A consequence is that a triangulation of any oriented surface can 
be encoded by cutting the surface along non-contractible cycles and see the surface as a 
planar map with boundaries. This method is a first attempt to generalize PS algorithm 
to higher genus. The obtained algorithm is asymptotically optimal (in terms of number 
of bits) but it is not linear, nor bijective. 

The goal of this paper is to present a new generalization of PS algorithm to higher 
genus based on some strong structural properties. Applied on a well chosen orientation 
of a toroidal triangulation, what remains after the execution of the algorithm is a rooted 
unicellular map (which corresponds to the natural generalization of trees when going 
to higher genus) that can be encoded optimally using 3, 2451 n bits. Moreover, the 
algorithm can be performed in linear time and leads to a new bijection between toroidal 
triangulations and a particular type of unicellular maps. 

The two main ingredients that make PS algorithm work in an orientation of a planar 
map are minimality and accessibility of the orientation. Minimality means that there 
is no clockwise cycle. Accessibility means that there exists a root vertex such that all 
the vertices have an oriented path directed toward the root vertex. Given a : y —>■ N, 
an orientation of G is an a-orientation if for every vertex u G P its outdegree dA{v) 
equals a{v). The existence and uniqueness of minimal orientations in the plane is given 
by the following result of Felsner [ 7 ] (related to older results of Propp m and Ossona 
de Mendez dl) : the set of a-orientations of a given planar map carries a structure of 
distributive lattice. This gives the existence and uniqueness of a minimal a-orientation 
as soon as an a-orientation exists. Felsner’s result enables several applications of PS 
method to other kind of planar maps, see 13 Eli- In all these cases the accessibility of 
the considered a-orientations is a consequence of the natural choice of a, like in Poulalhon 
and Schaeffer’s original work [TS] where any 3-orientation of the inner edges of a planar 
triangulation is accessible for any choice of root vertex on the outer face. (Note that the 
conventions may differs in the literature: the role of outgoing and incoming edges are 
sometimes exchanged and/or the role of clockwise and counterclockwise.) 

For higher genus, the minimality can be obtained by the following generalization of 
Felsner’s result. The second author, Knauer and the third author m showed that on 
any oriented surface the set of orientations of a given map having the same homology 
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carries a structure of distributive lattice. Note that a has been removed here since it is 
captured by the homology (see Sectionj^for a brief introduction to homology). Note also 
that this result is equivalent to an older result of Propp m where the lattice structure is 
described in the dual setting. Since this result is very general, there is hope to be able to 
further generalize PS method to surfaces. Note that a given map on an oriented surface 
can have several a-orientations (for the same given a) that are not homologous. So the 
set of a-orientations of a given map is now partitioned into distributive lattices contrarily 
to the planar case where there is only one lattice (and thus only one minimal element). 
In the case of toroidal triangulations we manage to face this problem and maintain a 
bijection by recent results on the structure of 3-orientations of toroidal triangulations. 
We identify a special lattice (and thus a special minimal orientation) using the notion of 
Schnyder woods generalized to the torus by the second and third author in [11] (further 
generalized in [T2|, see also |14) for a unified presentation). 

The main issue while trying to extend PS algorithm to higher genus is the accessibil¬ 
ity. Accessibility toward the outer face is given almost for free in the planar case because 
of Euler’s formula that sums to a strictly positive value. For an oriented surface of genus 
g > 1 new difficulties occur. Already in genus 1 (the torus), even if the orientation 
is minimal and accessible PS algorithm can visit all the vertices but not all the angles 
of the map because of the existence of non-contractible cycles. We can show that the 
special minimal orientation that we choose has the nice property that this problem never 
occurs. In genus g > 2 things get even more difficult with separating non-contractible 
cycles that make having accessibility of the vertices already difficult to obtain. 

Another problem is to recover the original map after the execution of the algorithm. 
If what remains after the execution of PS method is an unicellular map then the map 
can be recovered with the same simple rules as in the plane. Unfortunately for many 
minimal orientations the algorithm leads to an unicellular embedded graph covering all 
the edges but that is not a map (the only face is not a disk) and it is not possible 
to directly recover the original map. Here again, the choice of our special orientation 
ensures that this never happens. 

Finally the method presented here can be implemented in linear time. Clearly the 
execution of PS algorithm is linear but the difficulty lies in providing the algorithm with 
the appropriate orientation in input. Computing the minimal Schnyder wood of a planar 
triangulation can be done in linear time quite easily by using a so-called shelling order 
(or canonical order, see m) Other similar ad-hoc linear algorithms can be found for 
other kinds of a-orientations of planar maps (see for example [9] Chapter 3]). Such 
methods are not known in higher genus. We solve this problems by first computing an 
orientation in our special lattice and then go down in the lattice to find the minimal 
orientation. All this can be performed in linear time. 

Generalizing the method presented here to higher genus and other kind of maps 
thus raises several challenging questions and we hope that the present paper will lead 
to further generalizations of planar bijections, coding, counting and sampling, to higher 
genus. 
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A brief introduction to homology and to the corresponding terminology used in the 
paper is given in Section [2j In Section we present the definitions and results we 
need concerning the generalization of Schnyder woods to the toroidal case. In Section]^ 
we introduce a reformulation of Poulalhon and Schaeffer’s original algorithm that is 
applicable to any orientation of any map on an oriented surface. The main theorem 
of this paper is proved in Section that is, for a toroidal triangulation given with an 
appropriate root and orientation, the output of the algorithm is a toroidal unicellular map 
covering all the edges of the triangulation. In Sectionj^ we show how one can recover the 
original triangulation from the output. This output is then used in Sectionj^to optimally 
encode a toroidal triangulation. The linear time complexity of the method is discussed in 
Section]^ In Section(resp. Section 11), we exhibit a bijection between appropriately 
rooted toroidal triangulations and rooted (resp. non-rooted) toroidal unicellular maps. 
To obtain the non-rooted bijection, further structural results concerning the particular 
Schnyder woods considered in this paper are given in Section 10 Finally, a possible 


generalization to higher genus is discussed in Section 12 


2 A bit of homology 

We need a bit of surface homology of general maps, which we discuss now. The pre¬ 
sentation is not standard but it is short and sufficient to fit our needs. For a deeper 
introduction to homology we refer to |in) . 

Consider a map G with edge set E, on an orientable surface of genus g, given with 
an arbitrary orientation of its edges. This fixed arbitrary orientation is implicit and is 
used to manipulate flows. A flow </> on G is a vector in Z®. For any e £ E, we denote 
by 4>e the coordinate e of cj). 

A walk VF of G is a sequence of edges such that consecutive edges are incident. A 
walk is closed if the starting and ending vertices are the same. A walk has a characteristic 
flow 4>{W) defined by: 

(j)(W)e ■= times W traverses e forward — times W traverses e backward 

This definition naturally extends to sets of walks. From now on we consider that a 
set of walks and its characteristic flow are the same object. We do similarly for oriented 
subgraphs as they can be seen as sets of walks. 

A facial walk is a closed walk bounding a face. Let E be the set of counterclockwise 
facial walks and let F =< 4i{E) > the subgroup of Z-® generated by E. Two flows (j), 
are said to be homologous (p — £ ¥. A flow (p is 0-homologous if it is homologous to 

the zero flow, i.e. (p £¥. 

Let W be the set of closed walks and let W =< ^(W) > the subgroup of Z® 
generated by W. The group H{G) = W/F is the first homology group of G. Since 
dim(W) = m — n -|- 1 and dim(F) = / — 1, Euler’s Formula gives dim(LI(G)) = 2g. So 
H{G) = 1?^ only depends on the genus of the map. A set (i?i,..., B 2 g) of (closed) walks 
of G is said to be a basis for the homology if {(p{Bi )^..., (p{B 2 g)) is a basis of H{G). 
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3 Toroidal Schnyder woods 


Schnyder [ 2 T] introduced Schnyder woods for planar triangulations using the following 
local property: 

Given a map G, a vertex v and an orientation and coloring of the edges incident to 
V with the colors 0, 1, 2, we say that a vertex v satishes the Schnyder property if (see 

Vertex v has out-degree one in each color. 

The edges eo(v), ei(v), 62 (u) leaving v in colors 0, 1, 2, respectively, occur in 
counterclockwise order. 

Each edge entering v in color i enters v in the counterclockwise sector from ej_|_i(u) 
to ei-i(v) (where i + 1 and i — 1 are understood modulo 3). 

1 



Figure 1: The Schnyder property. The correspondence between red, blue, green and 0, 
1 , 2 and the arrow shapes used here serves as a convention for all figures in the paper. 

Given a planar triangulation G, a (planar) Schnyder wood of G is an orientation and 
coloring of the inner edges of G with the colors 0, 1, 2 , where each inner vertex v satishes 
the Schnyder property. In mm (see also the HDR thesis of the third author |14] 1 a 
generalization of Schnyder woods for higher genus has been proposed. Since this paper 
deals with triangulations of the torus only, we use a simplihed version of the dehnitions 
and results from mm- 

The dehnition of Schnyder woods for toroidal triangulations is the following. Given a 
toroidal triangulation G, a (toroidal) Schnyder wood of G is an orientation and coloring 
of the edges of G with the colors 0 , 1 , 2 , where each vertex satishes the Schnyder property 
(see Figure]^ for an example). The three colors 0, 1, 2 are completely symmetric in the 
dehnition, thus we consider that two Schnyder woods that are obtained one from the 
other by a (cyclic) permutation of the colors are in fact the same object. We consider 
that a Schnyder wood and its underlying orientation are the same object since one can 
easily recover a coloring of the edges in a greedy way (by choosing the color of an edge 
arbitrarily and then satisfying the Schnyder property at every vertex). 
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Figure 2: A Schnyder wood of a toroidal triangulation (opposite sides are identified in 
order to form a torus). 


Note that the situation is quite different from the planar case. In a Schnyder wood 
of a toroidal triangulation, each vertex has exactly one outgoing arc in each color, so 
there are monochromatic cycles contrarily to the planar case (one can show that these 
monochromatic cycles are not contractible). Moreover the graph induced by one color is 
not necessarily connected. Also, by a result of De Fraysseix and Ossona de Mendez [8], 
there is a bijection between orientations of the internal edges of a planar triangulation 
where every inner vertex has outdegree 3 and Schnyder woods. Thus, in the planar case, 
any orientation with the proper outdegrees corresponds to a Schnyder wood. This is not 
true for toroidal triangulations since there exists 3-orientations that do not correspond 
to a Schnyder wood (see Figure]^. 



Figure 3: Two different orientations of a toroidal triangulation. Only the one on the 
right corresponds to a Schnyder wood. 

A Schnyder wood of a toroidal triangulation is said to be crossing, if for each pair i,j 
of different colors, there exists a monochromatic cycle of color i intersecting a monochro¬ 
matic cycle of color j. The existence of crossing Schnyder woods is proved in m Theo¬ 
rem 1] (note that in [11] the crossing property is included in the definition of Schnyder 
woods, see m for a unified presentation): 
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Theorem 1 ( jll| ) A toroidal triangulation admits a crossing Schnyder wood. 

Figure [^depicts two different Schnyder woods of the same graph where just the one 
on the left is crossing (on the right case the red and green monochromatic cycles do not 
intersect, we say that the Schnyder wood is “half-crossing” since blue crosses both green 
and red, see uadi! for a formal definition). Note that the Schnyder wood on the right is 
obtained from the one on the left by flipping a clockwise triangle into a counterclockwise 
triangle. 



Figure 4: A crossing and an half-crossing Schnyder wood. 


Consider a toroidal triangulation G given with a crossing Schnyder wood. Let Dq be 
the corresponding 3-orientation of G. Let 0{G) be the set of all the orientations of G 
that are homologous to L^o- A consequence of m Theorem 5 and Corollary 2] is that all 
the crossing Schnyder woods of G are homologous to each other. So 0{G) contains all 
the crossing Schnyder woods of G. Thus the definition of 0{G) does not depend on the 
particular choice of Dq and thus it is uniquely defined. A consequence of |12( Theorem 4 
and Corollary 2] is that every orientation of 0{G) corresponds to a Schnyder wood. 
Thus we call the elements of 0{G) the homologous-to-crossing Schnyder woods (or HTC 
Schnyder woods for short). Note that all the crossing Schnyder woods are HTC. 

Figurej^gives an example of an HTC Schnyder wood that is not crossing and a Schny¬ 
der woods that is not HTC. The example on the left is obtained from the crossing Schny¬ 
der wood of Figure]^ by flipping two triangles (one to obtain the half-crossing Schnyder 
wood of Figure and then another one flipped from counterclockwise to clockwise). 
Thus it is HTC since the difference with a crossing Schnyder wood is a 0-homologous 
oriented subgraph. The example on the right of Figure is obtained from the crossing 
Schnyder wood of Figure|^by reversing the three vertical red monochromatic cycles. The 
union of these three cycles is not a 0-homologous oriented subgraph, thus the resulting 
orientation is not HTC. 




Not HTC 


Figure 5: Non-crossing Schnyder woods. 
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It is proved in [12] that on any oriented surface the set of orientations of a given map 
having the same homology carries a structure of distributive lattice for a particular order 
defined below. Thus in particular the set of HTC Schnyder woods carries a structure of 
distributive lattice. 

Let us define an order on the orientations of G. For that purpose, choose an arbitrary 
face /o of G and let Fq be its counterclockwise facial walk (this choice of a particular 
face corresponds to the choice of the outer face in the planar case). Let F be the set 
of counterclockwise facial walks of G and F' = F \ Fq. We say that a 0-homologous 
oriented subgraph T of G is counterclockwise (resp. elockwise) w.r.t. /o, if its character¬ 
istic flow can be written as a combination with positive (resp. negative) coefficients of 
characteristic flows of F', i.e. <^(r) = YIfgT' with A G (resp. —A G NI-^'I). 

Given two orientations D and D' of G, let D\D' denote the subgraph of D induced 
by the edges that are not oriented as in D'. We set D <f^ D' if and only \{ D\D' is 
counterclockwise. In [121 Theorem 7] the following is proved: 

Theorem 2 ( {12| ) {0(G), <fg) is a distributive lattice. 

Since (0(G), </p) is a distributive lattice, it has a unique minimal element. The 
following lemma gives a property of this minimum that is essential to apply Poulalhon 
and Schaeffer’s method. 

Lemma 1 The minimal element of {0{G), <ff^) is the only HTC Schnyder wood that 
eontains no elockwise (non-empty) 0-homologous oriented subgraph w.r.t. /q. 

Proof. Let Omin be the minimal element of (0(G), </o)- Suppose by contradiction that 
Omin contains a clockwise non-empty 0-homologous oriented subgraph T w.r.t. /q. The 
orientation of G obtained from Omin by reversing all the edges of T gives an orientation 
D G 0(G) such that T = \ D. Furthermore, by definition of <jp, we have D <fQ 

Ojnin, a contradiction to the minimality of Z?min- So Dmin contains no clockwise non¬ 
empty 0-homologous oriented subgraph w.r.t. /q. 

We now show that this characterizes D mw . For any D G 0(G), distinct from D m\n . 
we have L>min </o D. Thus T = D \ iAmin is a non-empty clockwise 0-homologous 
oriented subgraph of iA. □ 

The crossing Schnyder wood of Figure]^ is the minimal HTC Schnyder wood for the 
choice of /o corresponding to the shaded face. This example is used in the next sections 
to illustrate Poulalhon and Schaeffer’s method. 

The two HTC Schnyder woods of Figure]^ are not minimal (for any choice of special 
face /o) since they contain several triangles that are oriented clockwise. On the contrary, 
the HTC Schnyder wood of Figure]^ is minimal w.r.t to its only face oriented clockwise. 
These examples shows that the minimal HTC Schnyder wood is not always crossing. 

We define the dual orientation D* of an orientation iA of G as an orientation of the 
edges of the dual map G* of G satisfying the following rule: the dual e* of an edge e 
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Figure 6: The minimal HTC Schnyder wood of Kj w.r.t. the shaded face. 

goes from the face on the left of e to the face on the right of e. The following lemma 
gives the key property of HTC Schnyder woods that we need in this paper: 

Lemma 2 If D is an orientation corresponding to an HTC Schnyder wood, then the 
dual orientation D* contains no oriented non-contractible cycle. 

Proof. We first prove the property for a crossing Schnyder wood and then show that 
it is stable by reversing a 0-homologous oriented subgraph. Thus it is true for all HTC 
Schnyder woods. 

Consider a crossing Schnyder wood of G by Theorem and let Dq be the corre¬ 
sponding orientation. For i G {0,1,2}, let Ci be a monochromatic cycle of color i. 
By m Theorem 7], in a crossing Schnyder wood, the monochromatic cycles are not 
contractible and any two monochromatic cycles of different colors are not homologous 
and intersecting. Thus for i G {0,1,2}, the two cycles Ci-i and Ci+i form a basis Bi 
for the homology. By the Schnyder property, cycle Ci-i is crossing Ct (maybe several 
time) from left to right. Thus the homology of any closed curve can be expressed in at 
least one of the basis Bi with only positive coefficients. 

Snppose now by contradiction that Dq contains an oriented non-contractible cycle 
C*. Let i in {0,1, 2}, snch that C* is homologous to Xi-iCi-i + Aj+iCj+i with Aj_i > 0 
and Aj+i > 0. Then Cj+i is crossing C* at least once from left to right, contradicting 
the fact that C* is an oriented cycle of Dq. So Dq contains no oriented non-contractible 
cycle. 

Consider now an HTC Schnyder wood of G and let D be the corresponding orien¬ 
tation. Since D and Dq are both element of 0{G) they are homologous to each other. 
Let T be the 0-homologous oriented snbgraph of D such that T = D \ Dq. Thus Dq is 
obtained from D by reversing the edges of T. 
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Suppose by contradiction that D* contains an oriented non-contractible cycle C*. 
The oriented subgraph T is 0-homologous thus it intersects C* exactly the same number 
of time from right to left than from left to right. Since C* is oriented T cannot intersect 
it from left to right. So T does not intersect C* at all. Thus reversing T to go from 
D to Dq does not affect C*. Thus C* is an oriented non-contractible cycle of Dq, a 
contradiction. □ 

For the non-HTC Schnyder wood of Figure one can see that there is an horizon¬ 
tal oriented non-contractible cycle in the dual, so it does not satisfy the conclusion of 
Lemma[^ Note that this property is not a characterization of being HTC. Figureis a 
Schnyder wood that is not HTC but satisfies the conclusion of Lemma (we leave the 
reader check that this Schnyder wood is not HTC, it will be easier after Section and 
the definition of 7 ). 



Figure 7: A Schnyder wood that is not HTC but contains no oriented non-contractible 
cycle in the dual. 


4 Poulalhon and Schaeffer’s algorithm on oriented surfaces 

In this section we introduce a reformulation of Poulalhon and Schaeffer’s original algo¬ 
rithm. This version is more general in order to be applicable to any orientation of any 
map on an oriented surface. The execution slightly differs from the original formula¬ 
tion, even on planar triangulations. In [18], the authors first delete some outer edges 
of the triangulation before executing the algorithm. We do not consider some edges to 
be special here since we want to apply the algorithm on any surface but the core of the 
algorithm is the same. We show general properties of the algorithm in this section before 
considering toroidal triangulations in the forthcoming sections. 

Algorithm PS 

Input : An oriented map G, a root vertex vq and a root edge eo incident to vg- 
Output : An embedded graph U with stems. 

1. Let V := Vq, e := eg, U := 0. 

2. Let v' be the extremity of e different from v. 

Case 1 : e is non-marked and entering v. Add e to U and let u := u'. 


10 















Case 2 : e is non-marked and leaving v. Add a stem to U incident to v and corre¬ 
sponding to e. 

Case 3 : e is already marked and entering v. Do nothing. 

Case 4 : e is already marked and leaving v. Let v := v'. 

3. Mark e. 

4. Let e be the next edge around v in counterclockwise order after the current e. 

5. While {v, e) ^ {vq, cq) go back to 2. 

6. Return U. 


We insist on the fact that the output of Algorithm PS is a graph embedded on 
the same surface as the input map but that this embedded graph is not necessarily a 
map (i.e some faces may not be homeomorphic to open disks). In the following section 
we show that in our specific case the output U is an unicellular map. 


Consider any oriented map G on an oriented surface given with a root vertex vq and 
a root edge cq incident to vq. When Algorithm PS is considering a couple (u,e) we 
see this like it is considering the angle at v that is just before e in clockwise order. The 
particular choice of vq and eo is thus in fact a particular choice of a root angle oq that 
automatically defines a root vertex vq, a root edge cq, as well as a root face /q. From 
now on we consider that the input of Algorithm PS is an oriented map plus a root 
angle (without specifying the root vertex, face and edge). 


The angle graph of G, is the graph defined on the angles of G and where two angles 
are adjacent if and only if they are consecutive around a vertex or around a face. An 
execution of Algorithm PS can be seen as a walk in the angle graph. Figurej^illustrates 
the behavior of the algorithm corresponding to Case 1 to 4. In each case, the algorithm is 
considering the angle in top left position and depending on the marking of the edge and 
its orientation the next angle that is considered is the one that is the end of the magenta 
arc of the angle graph. The cyan edge of Case 1 represents the edge that is added to 
U by the algorithm. The stems of U added in Case 2 are not represented in cyan, in 
fact we will represent them later by an edge in the dual. Indeed seeing the execution of 
Algorithm PS as a walk in the angle graph enables us to show that Algorithm PS 
behaves exactly the same in the primal or in the dual map (as explained later). 

On Figure]^ we give an example of an execution of Algorithm PS on the orientation 
corresponding to the minimal HTC Schnyder wood of Kj of Figure 


Let o be a particular angle of the map G. It is adjacent to four other angles in 
the angle graph (see Figure 10). Let v,f be such that a is an angle of vertex v and 
face /. The next-vertex (resp. previous-vertex) angle of a is the angle appearing just 
after (resp. before) a in counterclockwise order around v. Similarly, the next-face (resp. 
previous-face) angle of a is the angle appearing just after (resp. before) a in clockwise 
order around /. These definitions enable one to orient consistently the edges of the 
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il non-marked 


Case 1 



non-marked 


Case 2 




il marked 


Y marked 


Case 3 Case 4 

Figure 8: The four cases of Algorithm PS 




Figure 9: An execution of Algorithm PS on Kj given with the orientation correspond¬ 
ing to the minimal HTC Schnyder wood of Figure Vertices are numbered in black. 
The root angle is identified by a root symbol and chosen in the face for which the orien¬ 
tation is minimal (i.e. the shaded face of Figure]^. The magenta arrows and numbers 
are here to help the reader to follow the cycle in the angle graph. The output U is a 
toroidal unicellular map, represented here as an hexagon where the opposite sides are 
identihed. 


angle graph like in Figure 10 so that for every oriented edge {a, a'), a' is a next-vertex 
or next-face angle of a. 

The different cases depicted in Figure show that an execution of Algorithm PS 
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Figure 10: Orientation of the edges of the angle graph. 


is just an oriented walk in the angle graph (i.e. a walk that is following the orientation 
of the edges described in Figure [T^. The condition in the while loop ensures that when 
the algorithm terminates, this walk is back to the root angle. The following proposition 
shows that the algorithm actually terminates: 

Proposition 1 Consider an oriented map G on an oriented surface and a root angle 
oq. The execution o/A lgorithm PS on (G, uq) terminates and corresponds to a cycle 
in the angle graph. 

Proof. We consider the oriented walk W in the angle graph corresponding to the 
execution of Algorithm PS . Note that W may be infinite. The walk W starts with 
ao, and if it is hnite it ends with oq and contains no other occurrence of ao (otherwise the 
algorithm should have stopped earlier). Toward a contradiction, suppose that W is not 
simple (i.e. some angles different from the root angle ao are repeated). Let a 7 ^ ao be the 
first angle along W that is met for the second time. Let ai,a 2 be the angles appearing 
before the first and second occurrence of a in W, respectively. Note that ai / 02 by the 
choice of a. 

If ai is the previous-vertex angle of a, then 02 is the previous-face angle of a. When 
the algorithm considers ai, none of a and 02 are already visited, thus edge e is not 
marked. Since the execution then goes to a after ai, we are in Case 2 and the edge e 
between a and ai is oriented from v, where v is the vertex incident to a. Afterward, 
when the algorithm reaches 02 , Case 3 applies and the algorithm cannot go to a, a 
contradiction. The case where ai is the previous-face angle of a is similar. 

So W is simple. Since the angle graph is finite, W is finite. So the algorithm 
terminates, thus W ends on the root angle and kF is a cycle. □ 

In the next section we see that in some particular cases the cycle in the angle graph 
corresponding to the execution of PS algorithm (Proposition can be shown to be 
Hamiltonian like on Figure 

By Proposition an angle is considered at most once by Algorithm PS. This 
implies that the angles around an edge can be visited in different ways depicted on 
Figure pT| Consider an execution of Algorithm PS on G. Let C be the cycle formed 
in the angle graph by Proposition Let P be the set of edges of the output U (without 
the stems) and Q be the set of dual edges of edges of G corresponding to stems of U. 
These edges are represented on Figure [TT] in cyan for P and in yellow for Q. They are 
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considered with their orientation (recall that the dual edge e* of an edge e goes from the 
face on the left of e to the face on the right of e). Note that C does not cross an edge 
of P or Q, and moreover P and Q do not intersect (i.e. an edge can be in P or its dual 
in Q but both cases cannot happen). 



Case 1 alone Case 1 and then 4 




Figure 11: The different cases of Algorithm PS seen in a dual way. The number of 
the angles gives the order in which the algorithm visits them (unvisited angles are not 
numbered). The edges of P and Q are respectively cyan and yellow. 


One can remark that the cases of Figure im are dual of each other. One can see 
that Algorithm PS behaves exactly the same if applied on the primal map or on the 
dual map. The only modifications to make is to start the algorithm with the face /o as 
the root vertex, the dual of edge cq as the root edge and to replace counterclockwise by 
clockwise at Line 4. Then the cycle C formed in the angle graph is exactly the same and 
the output is Q with stems corresponding to P (instead of P with stems corresponding 
to Q). Note that this duality is also illustrated by the fact that the minimality of the 
orientation of G w.r.t. the root face is nothing else than the accessibility of the dual 
orientation toward the root face. Indeed, a clockwise 0-homologous oriented subgraph 
of G w.r.t /o corresponds to a directed cut of the dual where all the edges are oriented 
from the part containing /q. The following lemma shows the connectivity of P and Q: 


Lemma 3 At each step of the algorithm, for every vertex v appearing in an edge of P 
(resp. Q), there is an oriented path from v to vq (resp. fo) consisting only of edges of 
P (resp. Q). In particular P and Q are connected. 

Proof. If at a step a new vertex is reached then it correspond to Case 1 and the 
corresponding edge is added in P and oriented from the new vertex, so the property is 
satisfied by induction. As observed earlier the algorithm behaves similarly in the dual 
map. □ 
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Let C be the set of angles of G that are not in C. Any edge of G is bounded by 
exactly 4 angles. Since C is a cycle, the 4 angles around an edge are either all in C, all 
in C or 2 in each set (see Figure 11). Moreover, if they are 2 in each set, these sets are 
separated by an edge of P or an edge of Q. Hence the frontier between G and C is a 
set of edges of P and Q. Moreover this frontier is an union of oriented closed walks of 
P and of oriented closed walks of Q. In the next section we study this frontier in more 
details to show that C is empty in the case considered there. 


5 From toroidal triangulations to unicellular maps 

Let G be a toroidal triangulation. In order to choose appropriately the root angle oq, 
we have to consider separating triangles. A triangle is a closed walk of size 3 (it is 
not necessarily a cycle since non-contractible loops are allowed and it is not necessarily 
contractible). A separating triangle is a contractible triangle that is different from a face 
of G. We say that an angle is in the strict interior of a separating triangle if it is in its 
contractible region and not incident to a vertex of the triangle. We choose as root angle 
ao any angle that is not in the strict interior of a separating triangle. One can easily see 
that such an angle ao always exists. Indeed the interiors of two contractible triangles 
are either disjoint or one is included in the other. So, the angles that are incident to a 
contractible triangle whose interior is maximal by inclusion satisfy the property. 

The goal of this section is to prove the following theorem (see Figure]^ for an exam¬ 
ple): 

Theorem 3 Consider a toroidal triangulation G, a root angle ao that is not in the strict 
interior of a separating triangle and the orientation of the edges of G corresponding to 
the minimal HTC Schnyder wood w.r.t. the root face fo containing ao- Then the output 
U o/A lgorithm PS applied on (G, oq) is a toroidal unicellular map covering all the 
edges of G. 

Consider a toroidal triangulation G, a root angle ao that is not in the strict interior of 
a separating triangle and the orientation of the edges of G corresponding to the minimal 
HTC Schnyder wood w.r.t. the root face fo containing ao- Let P be the output of 
Algorithm PS applied on (G,ao). We use the same notation as in previous section: 
the cycle in the angle graph is G, the set of angles that are not in G is G, the set of 
edges of P is P, the dual edges of stems of P is Q- 

Lemma 4 The frontier between C and C contains no oriented closed walk of Q- 


Proof- Suppose by contradiction that there exists such a walk W- Then along this 
walk, all the dual edges of W are edges of G oriented from the region containing G 
toward G as one can see in Figure 
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If W is non-contractible, then W contains an 
oriented non-contractible cycle, a contradiction to Lemma So W is contractible. So 
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it contains an oriented contractible cycle W^ and then either C is in the contractible 
region delimited by W'^ or not. The two case are considered below: 

Suppose first that C lies in the non-contractible region of W. Then consider the 
plane map G' obtained from G by keeping only the vertices and edges that lie (strictly) 
in the contractible region delimited by W. Let n' be the number of vertices of G'. All 
the edges incident to G' that are not in G' are entering G'. So in G' all the vertices have 
outdegree 3 as we are considering 3-orientations of G. Thus the number of edges of G' 
is exactly 3n', contradicting the fact that the maximal number of edges of planar map 
on n vertices is 3n — 6 by Euler’s formula. 

Suppose now that C lies in the contractible region of W. All the dual edges of W 
are edges of G oriented from its contractible region toward its exterior. Consider the 
graph Gout obtained from G by removing all the edges that are cut by W and all the 
vertices and edges that lie in the contractible region of W. As G is a map, the face 
of Gout containing W' is homeomorphic to an open disk. Let F be its facial walk (in 
Gout) and let k be the length of F. We consider the map obtained from the facial walk 
F by putting back the vertices and edges that lied inside. We transform this map into 
a plane map G' by duplicating the vertices and edges appearing several times in F, in 
order to obtain a triangulation of a cycle of length k. Let be the number of 

vertices, edges and faces of G'. Every inner vertex of G' has outdegree 3, there are no 
other inner edges, so the total number of edges of G' is m' = 3{n' — k) + k. All the 
inner faces have size 3 and the outer face has size k, so 2m' = 3(/' — 1) + k. By Euler’s 
formula n' — m! + f = 2. Combining the three equalities gives k = 3 and F is hence a 
separating triangle of G. This contradicts the choice of the root angle, as it should not 
lie in the strict interior of a separating triangle. □ 

A subgraph of a graph is spanning if it is covering all the vertices. An Hamiltonian 
cycle is a spanning cycle. 


Lemma 5 The cycle G is an Hamiltonian cycle of the angle graph, all the edges of G 
are marked exactly twice, the subgraph Q of G* is spanning, and, ifn>2, the subgraph 
P of G is spanning. 


Proof. Suppose for a contradiction that C is non empty. By Lemma and Section 
the frontier T between G and G is an union of oriented closed walks of P. Hence a 
face of G has either all its angles in C or all its angles in G. Moreover T is a non¬ 
empty union of oriented closed walk of P that are oriented clockwise according to the 
set of faces containing G (see the hrst case of Eigure 11). This set does not contain 
/o since uq is in /o and G. As in Section let F be the set of counterclockwise facial 
walks of G and Fq be the counterclockwise facial walk of /q. Let F' = F \ Fq, and 
F^ C F' be the set of counterclockwise facial walks of the faces containing G. We have 
<f{T) = — ’YliF&T-^^iF). So T is a clockwise non-empty 0-homologous oriented subgraph 
w.r.t. /q. This contradicts Lemma[^and the minimality of the orientation w.r.t. /q. So 
G is empty, thus G is Hamiltonian and all the edges of G are marked twice. 
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Suppose for a contradiction that n > 2 and P is not spanning. Since the algorithm 
starts at uq, P is not covering a vertex u of G different from vq. Then the angles around 
V cannot be visited since by Figure [TT] the only way to move from an angle of one vertex 
to an angle of another vertex is through an edge of P incident to them. So P is spanning. 
The proof is similar for Q (note that in this case we have / > 2). □ 


Lemma 6 The first cycle created in P (resp. in Q) by the algorithm is oriented. 

Proof. Let e be the first edge creating a cycle in P while executing Algorithm PS 
and consider the steps of Algorithm PS before e is added to P. So P is a tree during 
all these steps. For every vertex of P we define P{v) the unique path from v to vq in 
P (while P is empty at the beginning of the execution, we define P{vq) = {uq}). By 
Lemma this path P{v) is an oriented path. We prove the following 

Claim 1 Consider a step of the algorithm before e is added to P and where the algorithm 
is considering a vertex v. Then all the angles around the vertices of P different from the 
vertices of P{v) are already visited. 

Proof. Suppose by contradiction that there is such a step of the algorithm where some 
angles around the vertices of P different from the vertices of P{v) have not been visited. 
Consider the first such step. Then clearly we are not at the beginning of the algorithm 
since P = P{v) = {uq}. So at the step just before, the conclusion holds and now it does 
not hold anymore. Clearly at the step before we were considering a vertex v' distinct 
from V, otherwise P{v) and P have not changed and we have the conclusion. So from v' 
to V we are either in Case 1 or Case 4 of Algorithm PS. If v has been considered by 
Case 1, then P{v) contains P{v') and the conclusion holds. If v has been considered by 
Case 4, then since P is a tree, all the angles around v' have been considered and v' is 
the only element of P \ P{v) that is not in P \ P{v'). Thus the conclusion also holds. □ 

Consider the iteration of Algorithm PS where e is added to P. The edge e is added 
to P by Case 1, so e is oriented from a vertex u to a vertex v such that v is already in 
P or u is the root vertex vq. Consider the step of the algorithm just before u is added 
to P. By Claim vertex u is not in P \ P{v) (otherwise e would have been considered 
before and it would be a stem). So u G Piv) and P{v) U {e} induces an oriented cycle 
of G. The proof is similar for Q. □ 


Lemma 7 P is a spanning unicellular map of G and Q is a spanning tree of G*. More¬ 
over one is the dual of the complement of the other. 

Proof. Suppose that Q contains a cycle, then by Lemma it contains an oriented 
cycle of G*. This cycle is contractible by Lemma Recall that by Lemma C is an 
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Hamiltonian cycle, moreover it does not cross Q, a contradiction. So Q contains no cycle 
and is a tree. 

By Lemma all the edges of G are marked at the end. So every edge of G is either 
in P or its dual in Q (and not both). Thus P and Q are the dual of the complement of 
each other. So P is the dual of the complement of a spanning tree of G*. Thus P is a 
spanning unicellular map of G. □ 


Theorem]^ is then a direct reformulation of Lemmaby the definition of P and Q: 

A toroidal unicellular map on n vertices has exactly n + 1 edges: n — 1 edges of a 
tree plus 2 edges corresponding to the size of a basis of the homology (i.e. plus 2g in 
general for an oriented surface of genus g). Thus a consequence of Theorem]^ is that the 
obtained unicellular map U has exactly n vertices, n +1 edges and 2n — 1 stems since the 
total number of edges is 3n. The orientation of G induces an orientation of U such that 
the stems are all outgoing, and such that while walking clockwise around the unique 
face of U from oq, the first time an edge is met, it is oriented counterclockwise according 
to this face, see Figure 12 where all the tree-like parts and stems are not represented. 
There are two types of toroidal unicellular maps depicted on Figure Two cycles of 


U may intersects either on a single vertex (square case) or on a path (hexagonal case). 
The square can be seen as a particular case of the hexagon where one side has length 
zero and thus the two corners of the hexagon are identihed. 




Figure 12: The two types of rooted toroidal unicellular maps. 


On Figure [T^ we give several examples of executions of Algorithm PS on minimal 
3-orientations. These examples show how important is the choice of the minimal HTC 
Schnyder wood in order to obtain Theorem In particular, the third example shows 
that Algorithm PS can visit all the angles of the triangulation (i.e. the cycle in the 
angle graph is Hamiltonian) without outputting an unicellular map. 

Note that the orientations of Figure are not Schnyder woods. One may wonder if 
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Figure 13: Examples of minimal 3-orientations that are not HTC Schnyder woods and 
where Algorithm PS respectively: (1) does not visit all the vertices, (2) visits all 
the vertices but not all the angles, and (3) visits all the angles but does not output an 
unicellular map. 


the fact of being a Schnyder wood is of any help for our method. This is not the case since 
there are examples of minimal Schnyder woods that are not HTC and where Algorithm 
PS does not visit all the vertices. One can obtain such an example by replicating 3 
times horizontally and then 3 times vertically the second example of Figure to form 
a 3 X 3 tiling and starts Algorithm PS from the same root angle. Conversely, there 
are minimal Schnyder woods that are not HTC where Algorithm PS does output a 
toroidal unicellular map covering all the edges (the Schnyder wood of Figure]^ can serve 
as an example while starting from an angle of the only face oriented clockwise). 


6 Recovering the original triangulation 

This section is dedicated to show how to recover the original triangulation from the 
output of Algorithm PS. The method is very similar to [18| since like in the plane the 
output has only one face that is homeomorphic to an open disk (i.e. a tree in the plane 
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and an unicellular map in general). 


Theorem 4 Consider a toroidal triangulation G, a root angle oq that is not in the strict 
interior of a separating triangle and the orientation of the edges of G corresponding to 
the minimal HTC Schnyder wood w.r.t. the root face /o containing qq- From the output 
U o/ Algorithm PS applied on {G,ao) one can reattach all the stems to obtain G by 
starting from the root angle oq and walking along the face of U in counterclockwise order 
(according to this face): each time a stem is met, it is reattached in order to create a 
triangular face on its left side. 


Theorem is illustrated on Figure where one can check that the obtained toroidal 
triangulation is K-/ (like on the input of Figure]^. 




output of Algorithm PS 


In fact in this section we define a method, more general than the one described in 
Theorem]^ that will be useful in next sections. 

Let Ur{n) denote the set of toroidal unicellular maps U rooted on a particular angle 
ao, with exactly n vertices, n + 1 edges and 2n — 1 stems satisfying the following: every 
vertex has exactly 2 stems, except the root vertex vq that has 1 more stem, and if the 
map is hexagonal, the two corners that have 1 less stem each, and if the map is square, 
the only corner that has 2 less stems (if the root vertex is a corner we simply combine 
the more and less). Note that the output U of Algorithm PS given by Theorem]^ is 
an element of Ur (n). 

Similarly to the planar case [18], we define a general way to reattached step by step 
all the stems of an element U of Ur{n). Let Uq = U, and, for 1 < A: < 2n — 1, let 
be the map obtained from Uk-i by reattaching one of its stem (we explicit below which 
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stem is reattached and how). The special face of Uq is its only face. For 1 < k < 2n — 1, 
the special face of Uk is the face on the right of the stem of Uk-i that is reattached 
to obtain Uk- For 0 < k < 2n — 1, the border of the special face of Uk consists of a 
sequence of edges and stems. We define an admissible triple as a sequence (ei,e 2 ,s), 
appearing in counterclockwise order along the border of the special face of Uk , such that 
ei = {u,v) and 62 = {v,w) are edges of Uk and s is a stem attached to w. The closure 
of the admissible triple consists in attaching s to u, so that it creates an edge {w, u) 
oriented from w to u and so that it creates a triangular face {u, v, w) on its left side. 
The complete closure of U consists in closing a sequence of admissible triple, i.e. for 
1 < k < 2n — 1, the map Uk is obtained from Uk-i by closing any admissible triple. 

Note that, for 0 < A: < 2n — 1, the special face of Uk contains all the stems of Uk- 
The closure of a stem reduces the number of edges on the border of the special face and 
the number of stems by 1 . At the beginning, the unicellular map Uq has n + 1 edges 
and 2n — 1 stems. So along the border of its special face, there are 2 n + 2 edges and 
2n — 1 stems. Thus there is exactly three more edges than stems on the border of the 
special face of Uq and this is preserved while closing stems. So at each step there is 
necessarily at least one admissible triple and the sequence Uk is well defined. Since the 
difference of three is preserved, the special face of U 2 n -2 is a quadrangle with exactly 
one stem. So the reattachment of the last stem creates two faces that are triangles and 
at the end U 2 n-i is a toroidal triangulation. Note that at a given step there might be 
several admissible triples but their closure are independent and the order in which they 
are performed does not modify the obtained triangulation U 2 n-i- 

We now apply the closure method to our particular case. Consider a toroidal trian¬ 
gulation G, a root angle oq that is not in the strict interior of a separating triangle and 
the orientation of the edges of G corresponding to the minimal HTC Schnyder wood 
w.r.t. the root face /q. Let U be the output of Algorithm PS applied on (( 7 , 00 ). 

Lemma 8 When a stem of U is reattached to form the corresponding edge of G, it splits 
the (only) face ofU into two faces- The root angle ofU is in the face that is on the right 
side of the stem- 

Proof- By Lemma the execution of Algorithm PS corresponds to an Hamiltonian 
cycle G = {oq, ... ,a 2 m;“o) in the angle graph of G- Thus G defines a total order < 
on the angles of G where ai < aj if and only if f < j- Let us consider now the angles 
on the face of U- Note that such an angle corresponds to several angles of G, that 
are consecutive in C and that are separated by a set of incoming edges of G (those 
incoming edges corresponding to stems of U)- Thus the order on the angles of G defines 
automatically an order on the angles of U- The angles of U considered in clockwise order 
along the border of its face, starting from the root angle, correspond to a sequence of 
strictly increasing angles for <. 

Consider a stem s of U that is reattached to form an edge e of G. Let Og be the 
angle of U that is situated just before s (in clockwise order along the border of the face 
of U) and a( be the angle of U where s should be reattached. If < a^, then when 
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Algorithm PS consider the angle a*, the edge corresponding to s is already marked 
and we are not in Case 2 of Algorithm PS . So < a'g and ao is on the right side of 
s. □ 

Recall that U is an element of Ur{n) so we can apply on U the complete closure 
procedure described above. We use the same notation as before, i.e. let Uq = U and for 
1 < k < 2n — 1, the map Uk is obtained from U^-i by closing any admissible triple. The 
following lemma shows that the triangulation obtained by this method is G: 

Lemma 9 The complete closure of U is G, i.e. U 2 n-i = G. 

Proof. We prove by induction on k that every face of Uk is a face of G, except for the 
special face. This is true for k = 0 since Uq = U has only one face, the special face. Let 
0 < A; < 2n — 2, and suppose by induction that every non-special face of Uk is a face of 
G. Let (ei, 62 , s) be the admissible triple of Uk such that its closure leads to Uk+i, with 
ei = {u,v) and 62 = iv,w). The closure of this triple leads to a triangular face {u,v,w) 
of Uk+i. This face is the only “new” non-special face while going from Uk to Uk+i. 

Suppose, by contradiction, that this face (n, v, w) is not a face of G. Let a„ (resp. 
ttw) be the angle of Uk at the special face, between ei and 62 (resp. 62 and s). Since G 
is a triangulation, and (tt, v, w) is not a face of G, there exists at least one stem of Uk 
that should be attached to Oy or Ow to form a proper edge of G. Let s' be such a stem 
that is the nearest from s. In G the edges corresponding so s and s' should be incident 
to the same triangular face. Let x be the origin of the stem s'. Let 2 ; G such 

that s' should be reattached to z. If 2 ; = v, then s should be reattached to x to form a 
triangular face of G. If 2 ; = w, then s should be reattached to a common neighbor of w 
and X located on the border of the special face of Uk in counterclockwise order between 
w and X. So in both cases s should be reattached to a vertex y located on the border of 
the special face of Uk in counterclockwise order between w and x (with possibly y = x). 
To summarize s goes from w to y and s' from x to 2 ;, and 2 ;, x, y, w appear in clockwise 
order along the special face of Uk- By Lemmathe root angle is on the right side of 
both s and s', this is not possible since their right sides are disjoint, a contradiction. 

So for 0 < /c < 2n — 2, all the non-special faces of Uk are faces of G. In particular 
every face of U 2 n-i except one is a face of G. Then clearly the (triangular) special face 
of U 2 n-i is also a face of G, hence U 2 n-i = G. □ 

Lemma shows that one can recover the original triangulation from U with any 
sequence of admissible triples that are closed successively. This does not explain how 
to find the admissible triples efficiently. In fact the root angle can be used to find a 
particular admissible triple of Uk'. 

Lemma 10 ForO < k < 2n—2, let s be the first stem met while walking counterclockwise 
from oo in the special face of Uk. Then before s, at least two edges are met and the last 
two of these edges form an admissible triple with s. 
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Proof. Since s is the first stem met, there are only edges that are met before s. Suppose 
by contradiction that there is only zero or one edge met before s. Then the reattachment 
of s to form the corresponding edge of G is necessarily such that the root angle is on the 
left side of s, a contradiction to Lemma So at least two edges are met before s and 
the last two of these edges form an admissible triple with s. □ 


Lemma 10 shows that one can reattach all the stems by walking once along the face 
of U in counterclockwise order. Thus we obtain Theorem 01 


Note that U is such that the complete closure procedure described here never wraps 
over the root angle, i.e. when a stem is reattached, the root angle is always on its right 
side (see Lemmaj^. The property of never wrapping over the root angle is called balanced 
in [2]. Let Virgin) denote the set of elements oiUr{n) that are balanced. So the output U 
of Algorithm PS given by Theoremj^is an element oilAr^hif^)- We exhibit in Section]^ 
a bijection between appropriately rooted toroidal triangulations and a particular subset 
of Urfiin). 


The possibility to close admissible triples in any order to recover the original trian¬ 
gulation is interesting comparing to the simpler method of Theorem since it enables 
to recover the triangulation even if the root angle is not given. This property is used 
in Section 11 to obtain a bijection between toroidal triangulations and some unrooted 
unicellular maps. 


Moreover if the root angle is not given, then one can simply start from any angle 
of U, walk twice around the face of U in counterclockwise order and reattached all the 
admissible triples that are encountered along this walk. Walking twice ensure that at 
least one complete round is done from the root angle. Since only admissible triples are 
considered, we are sure that no unwanted reattachment is done during the process and 
that the final map is G. This enables to reconstruct G in linear time even if the root 
angle is not known. This property is used in Section 


7 Optimal encoding 

The results presented in the previous sections allow us to generalize the encoding of 
planar triangulations, defined by Poulalhon and Schaeffer jT8|, to triangulations of the 
torus. The construction is direct and it is hence really different from the one of [5] 
where triangulations of surfaces are cut in order to deal with planar triangulations with 
boundaries. Here we encode the unicellular map outputted by Algorithm PS by a 
plane rooted tree with n vertices and with exactly two stems attached to each vertex, plus 
0(log(n)) bits. As in [5], this encoding is asymptotically optimal and uses approximately 
3.25n bits. The advantage of our method is that it can be implemented in linear time. 
Moreover we believe that our encoding gives a better understanding of the structure of 
triangulations of the torus. It is illustrated with news bijections that are obtained in 
Sections M and El 

Consider a toroidal triangulation G, a root angle oq that is not in the strict interior of 
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a separating triangle and the orientation of the edges of G corresponding to the minimal 
HTC Schnyder wood w.r.t. the root face /o- Let U be the output of Algorithm 
PS applied on (G, oo). As already mentioned at the end of Section]^ to retrieve the 
triangulation G one just needs to know U without the information of its root angle (by 
walking twice around the face of U in counterclockwise order and reattached all the 
admissible triples that are encountered along this walk, one can recover G). Hence to 
encode G, one just has to encode U without the position of the root angle around the 
root vertex (see Figure 15,(a)). 

By Lemma the unicellular map U contains a spanning tree T which is oriented 
from the leaves to the root vertex. The tree T contains exactly n — 1 edges, so there is 
exactly 2 edges of U that are not in T. We call these edges the special edges of U. We cut 
these two special edges to transform them into stems of T (see Figures 15 (a) and (b)). 
We keep the information of where are the special stems in T and on which angle of T 
they should be reattached. This information can be stored with 0(log(n)) bits. One 
can recover U from T by reattaching the special stems in order to form non-contractible 
cycles with T (see Figure [15} (c)). 





Figure 15: From unicellular maps to trees with special stems and back. 


So T is a plane tree on n vertices, each vertex having 2 stems except the root vertex 
Vo having three stems. Choose any stem so of the root vertex, remove it and consider 
that T is rooted at the angle where sq should be attached. The information of the root 
enable to put back so at its place. So now we are left with a rooted plane tree T on n 
vertices where each vertex has exactly 2 stems (see Figure [Ib} (a)). 


This tree T can easily be encoded by a binary word on 6n — 2 bits: that is, walking 
in counterclockwise order around T from the root angle, writing a “1” when going down 
along T, and a “0” when going up along T (see Figure}^ (a)). As in [18], one can encode 
T more compactly by using the fact that each vertex has exactly two stems. Thus T 
is encoded by a binary word on 4n — 2 bits: that is, walking in counterclockwise order 
around T from the root angle, writing a “1” when going down along an edge of T, and 
a “0” when going up along an edge or along a stem of T (see Figure [l6|(b) where the 
“red I’s” of Figure 16 (a) have been removed). Indeed there is no need to encode when 
going down along stems, this information can be retrieved afterward. While reading the 
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binary word to recover T, when a “0” is met, we should go up in the tree, except if the 
vertex that we are considering does not have already its two stems, then in that case we 
should create a stem (i.e. add a “red 1” before the “0”). So we are left with a binary 
word on 4n — 2 bits with exactly n — 1 bits “1” and 3n — 1 bits “0”. 



1010111011101001010010110100010101101000 

(a) 



(b) 


Figure 16: Encoding a rooted tree with two stems at each vertex. 

Similarly to [T^, using [U Lemma 7], this word can then be encoded with a binary 
word of length log2 ~ ~ 3, 2451 n bits. Thus we have the 

following theorem whose linearity is discussed in Section 


Theorem 5 Any toroidal triangulation on n vertices, can he encoded with a binary word 
of length 3.2451n + o(n) bits, the encoding and decoding being linear in n. 


8 Linear complexity 

In this section we show that the encoding method described in this paper, that is encod¬ 
ing a toroidal triangulation via an unicellular map and recovering the original triangula¬ 
tion, can be performed in linear time. The only difficulty lies in providing Algorithm 
PS with the appropriate input it needs in order to apply Theorem Then clearly the 
execution of Algorithm PS, the encoding phase and the recovering of the triangulation 
are linear. Thus we have to show how one can hnd in linear time a root angle oq that is 
not in the strict interior of a separating triangle, as well as the minimal HTC Schnyder 
wood w.r.t. the root face /q. 

Consider a toroidal triangulation G. Let us see how one can build a Schnyder wood 
of G in linear time. The contraction of a non-loop-edge e of G is the operation consisting 
of continuously contracting e until merging its two ends, as shown on Figure HZl Note 
that only one edge of each pair of edges forming a contractible 2-cycle is preserved (edges 
Cwx and Cwy on the figure). 
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Figure 17; The contraction operation. 


An edge uv is said to be contractible if it is not a loop and if it is not on a separating 
triangle (i.e. if after contracting uv one obtains a triangulation that is still without 
contractible 1- or 2-cycles). In |llj the existence of crossing Schnyder wood is proved by 
contraction. Unfortunately this proof cannot easily be transformed into a linear algo¬ 
rithm because of the crossing property that has to be maintained during the contraction 
process. Nevertheless we use contractions to obtain non-necessarily crossing Schnyder 
woods. If the triangulation obtained after contracting a contractible edge admits a 
Schnyder wood it is then easy to obtain a Schnyder wood of G. The rules for decon¬ 
tracting an edge in the case of toroidal triangulations are depicted on m Figure 21] 
where for each case one can choose any of the proposed colorings. For any toroidal 
triangulation, one can find contractible edges until the toroidal map has only one vertex 
(see [E]). A Schnyder wood of the toroidal map on one vertex is depicted on the right 
of Figure Thus one can obtain a Schnyder wood of any toroidal triangulation by this 
process. Nevertheless, to maintain linearity we have to be more precise since it is not 
trivial to find contractible edges. 

Consider an edge uv of G with incident faces uvx and vuy such that these vertices 
appear in clockwise order around the corresponding face (so we are in the situation of 
Figure [TT]). If u and v have more common neighbors, then consider their second common 
neighbor going clockwise around u from uv (the first one being x, and the last being y) 
and call it x'. Call y' their second common neighbor going counterclockwise around u 
from uv. Then either uvx' or uvy' is a separating triangle or edge uv is contractible. We 
consider these two cases below: 

• If uv is contractible, then it is contracted and we apply the procedure recursively 
to obtain a Schnyder wood of the contracted graph. Then we update the Schnyder 
wood as described above. Note that this update is done in constant time. 

• If uvx' (resp. uvy') is a separating triangle, one can remove its interior, recursively 
obtain a toroidal Schnyder wood of the remaining toroidal triangulation, build a 
planar Schnyder wood of the planar triangulation inside uvx' (resp. uvy'), and then 
superimpose the two (by eventually permuting the colors) to obtain a Schnyder 
wood of the whole graph. Note that computing a planar Schnyder wood can be 
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done in linear time using a canonical ordering (see |13jl. 

The difficulty here is to test if uvx' or uvy' are contractible triangles. For that pur¬ 
pose, one first need to compute a basis (i?i, B 2 ) for the homology. Consider a spanning 
tree of the dual map G*. The map obtained from G by removing those edges is unicel¬ 
lular, and removing its treelike parts one obtains two cycles {Bi,B 2 ) (intersecting on a 
path with at least one vertex) that form a basis for the homology. This can be computed 
in linear time for G and then updated in constant time when some edge is contracted or 
when the interior of some separating triangle is removed. Then a closed walk VF, given 
with an arbitrary orientation, is contractible if and only if W crosses Bi from right to 
left as many times as W crosses Bi from left to right, for i G {1,2}. This test is linear 
in |VF| hence constant time for the triangles uvx' and uvy'. Vertex u is fixed during the 
whole process so the total running time to compute a Schnyder wood of G is linear. 

From this Schnyder wood, one can compute in linear time a root angle uq not in the 
strict interior of a separating triangle. First note that in a 3-orientation of a toroidal 
triangulation, the edges that are inside a separating triangle and that are incident to 
the three vertices on the border are all oriented toward these three vertices by Euler’s 
formula. Thus an oriented non-contractible cycle cannot enter in the interior of a sepa¬ 
rating triangle. Now follow any oriented monochromatic path of the Schnyder wood and 
stop the first time this path is back to a previously met vertex vq. The end of this path 
forms an oriented monochromatic cycle G containing vo- If C is contractible then Euler’s 
formula is violated in the contractible region. Thus G is an oriented non-contractible 
cycle and cannot contain some vertices that are in the interior of a separating triangle. 
So Vq is not in the interior of a separating triangle and we can choose as root angle oq 
any angle incident to vq. 

In |12[ Section 9] it is proved how one can transform any 3-orientation (hence a 
Schnyder wood) of a toroidal triangulation into an half-crossing (hence HTC) Schnyder 
wood. The method consists in computing a so called “middle-path” (a directed path 
where the next edge chosen is the one leaving in the “middle”) and reversing some 
non-contractible “middle-cycles”. Clearly the method is linear even if not explicitly 
mentioned in [12]. Let Do be the corresponding obtained orientation of G. 

It remains to compute the minimal HTC Schnyder wood w.r.t. the root face /q. 
There is a generic known method m (see also eqi p.23]) to compute in linear time a 
minimal a-orientation of a planar map as soon as an a-orientation is given. This method 
also works on oriented surfaces and can be applied to obtain the minimal HTC Schnyder 
wood in linear time. We explain the method briefly below. 

It is much simpler to compute the minimal orientation Umin homologous to Dq in 
a dual setting. The first observation to make is that two orientations Di,D 2 of G are 
homologous if and only if there dual orientations D\ , D 2 of G* are equivalent up to 
reversing some directed cuts. Furthermore Di D 2 if and only if can be obtained 
from D 2 by reversing directed cuts oriented from the part containing /q. Let us compute 
-^min which is the only orientation of G*, obtained from Dq by reversing directed cuts, 
and without any directed cut oriented from the part containing /q. For this, consider 
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the orientation Dq of G* = {F,E*) and compute the set X C of vertices of G* that 
have an oriented path toward /o- Then {X,F \ X) is a directed cut oriented from the 
part containing /q that one can reverse. Then update the set of vertices that can reach 
/o and go on until X = F. It is not difficult to see that this can be done in linear time. 
Thus we obtain the minimal HTC Schnyder wood w.r.t. /o in linear time. 


9 Bijection with rooted unicellular maps 


Given a toroidal triangulation G with a root angle ao, we have defined a unique associated 
orientation: the minimal HTC Schnyder wood w.r.t. the root face /q. Suppose that G is 
oriented according to the minimal HTC Schnyder wood. If oq is not in the strict interior 
of a separating triangle then Theorems and show that the execution of Algorithm 
PS on {G, ao) gives a toroidal unicellular map with stems from which one can recover the 
original triangulation. Thus there is a bijection between toroidal triangulations rooted 
from an appropriate angle and their image by Algorithm PS. The goal of this section 
is to describe this image. 

Recall from Section that the output of Algorithm PS on (G, ao) is an element of 
Urchin). One may hope that there is a bijection between toroidal triangulations rooted 
from an appropriate angle and Ur^b{n) since this is how it works in the planar case. 
Indeed, given a planar triangulation G, there is a unique orientation of G (the minimal 
Schnyder wood) on which Algorithm PS, performed from an outer angle, outputs a tree 
covering all the edges of the graph. In the toroidal case, things are more complicated since 
the behavior of Algorithm PS on minimal HTC Schnyder woods does not characterize 
such orientations. 


Figure 18 gives an example of two (non-homologous) orientations of the same trian¬ 
gulation that are both minimal w.r.t. the same root face. For these two orientations, 
the execution of Algorithm PS from the same root angle gives two different elements 
of Urfi{2) (from which the original triangulation can be recovered by the method of 
Theorem]^. Thus we have to exhibit a particular property of HTC Schnyder woods 
that can be used to characterize which particular subset of Ur,b{n) is in bijection with 
appropriately rooted toroidal triangulations. 




Figure 18: A graph that can be represented by two different unicellular maps. 


For that purpose we use the following dehnition of 7 introduced in |12j . Consider a 
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particular orientation of G. Let C be a cycle that is given with an arbitrary direction 
{C is not necessarily a directed cycle). Then 7 ( 6 ') is defined by: 

7 ( 6 ') = # edges leaving C on its right — # edges leaving C on its left 


By the Schnyder property, it is clear that in a toroidal Schnyder wood, a monochro¬ 
matic cycle C always satisfies 7 ( 6 ') = 0. Consider a crossing Schnyder wood of G and 
C'i,C 2 two monochromatic cycles of different colors. Thus we have 7 (C'i) = 7 (C' 2 ) = 0. 
By m Theorem 7], the two cycles C'i,C '2 are non-contractible and non-homologous, 
thus they form a basis for the homology. While returning a 0-homologous oriented sub¬ 
graph, the value of 7 on a given cycle does not change. Thus any HTC Schnyder wood 
also satisfies 7 (C'i) = 7 (^ 2 ) = 0. Moreover it is proved in |12l Lemma 18] that if a 
3-orientation of a toroidal triangulation satisfies 7 equals 0 for two cycles forming a 
basis for the homology, then 7 equals 0 for any non-contractible cycle. Thus any HTC 
Schnyder wood satisfies 7 equals 0 for any non-contractible cycle. We call this property 
the 7 o property. Note that, for a 3-orientation, it is sufficient to satisfy 7 equals 0 on 
any two cycles forming a basis for the homology to have the 70 property. 


Actually the 70 property characterizes the 3-orientations that are HTC Schnyder 
woods. Indeed a consequence of |12| Theorem 5 and Lemma 18] is that if two 3- 
orientations both satisfy the 70 property, then they are homologous to each other and 
thus HTC. Note that for the 3-orientation on the right of Figure 18, we have 7 equals 
±2 for the horizontal cycle and this explain why this orientation is not HTC (one can 
find similar arguments for previous examples of non-HTC Schnyder woods presented in 
this paper, see Figures]^ and [^. 


Let us translate this 70 property on Ur{n). Consider an element U of Ur{n) whose 
edges and stems are oriented w.r.t. the root angle as follows: the stems are all outgoing, 
and while walking clockwise around the unique face of U from oq, the first time an edge 
is met, it is oriented counterclockwise w.r.t. the face of U. Then one can compute 7 on 
the cycles of U (edges and stems count). We say that an unicellular map oiUr{n) satisfies 
the 7 o property if 7 equals zero on its (non-contractible) cycles. Let us call 7 /r,b, 7 o (?t-) 
the set of elements of Urf ,{ n ) satisfying the 70 property. So the output of Algorithm 
PS given by Theorem[^is an element oi Ur^h,')Q{n)■ 

Let Tr{n) be the set of toroidal triangulations on n vertices rooted at an angle that 
is not in the clockwise interior of a separating triangle. Then we have the following 
bijection: 


Theorem 6 There is a bijection between Tr{n) and Ur 

Proof. Consider the mapping g that associates to an element of Tr{n), the output of 
Algorithm PS executed on the minimal HTC Schnyder wood w.r.t. the root face. By 
the above discussion the image of g is in Urfi^^g^n) and g is injective since one can recover 
the original triangulation from its image by Theorem]^ 

Conversely, given an element U of 7/^,6 ,70 (^) with root angle uq, one can build a 
toroidal map G by the complete closure procedure described in Section]^ The number 
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of stems and edges of U implies that G is a triangulation. Recall that oq defines an 
orientation on the edges and stems of U. Consider the orientation D ol G induced 
by this orientation. Since U is balanced, the execution of Algorithm PS on (G,ao) 
corresponds to the cycle in the angle graph of U obtained by starting from the root angle 
and walking clockwise in the face of U. Thus the output of Algorithm PS executed on 
(G, oo) is U. It remains to show that G is appropriately rooted and that D corresponds 
to the minimal HTC Schnyder wood w.r.t. this root, then g will be surjective. 

First note that by definition ofUr{n), the orientation D is a 3-orientation. 

Suppose by contradiction that oq is in the strict interior of a separating triangle. 
Then, since we are considering a 3-orientation, by Euler’s formula, the edges in the inte¬ 
rior of this triangle and incident to its border are all entering the border. So Algorithm 
PS started from the strict interior cannot visit the vertices on the border of the triangle 
and outside. Thus the output of Algorithm PS is not a toroidal unicellular map, a 
contradiction. So oq is not in the strict interior of a separating triangle. 

The 7 o property of U implies that 7 equals zero on two cycles of U. Hence these two 
cycles considered in G also satisfy 7 equals 0 and form a basis for the homology. So D 
is an HTC Schnyder wood. 

Suppose by contradiction that D is not minimal. Then, by Lemma it contains a 
clockwise ^on-empty) 0-homologous oriented subgraph w.r.t. /q. With the notations 
of Section let T be such a subgraph with 4'{T) = — YIf&F' with A G L 

Let Xfq = 0, and Amax = Xp. For 0 < i < A m a x. let Aj = {F G F\Xf > i}- 

For 1 < z < A max . let Ti be the oriented subgraph such that 0(ri) = —Y^F^Xi 
Then we have 4>{T) = Since T is an oriented subgraph, we have 

4>{T) G {—1,0, Thus for any edge of G, incident to faces Fi and F 2 , we have 

(Apj — XF 2 ) G {—1,0,1}. So, for 1 < z < A max . the oriented graph Ti is the frontier 
between the faces with A value equal to z and z — 1. So all the Tj are edge disjoint 
and are oriented subgraphs of D. Since T is non-empty, we have Amax > I, and Ti is 
non-empty. All the edges of Ti have a face of Xi on their right and a face of Xq on 
their left. Since U is an unicellular map, and Ti is a (non-empty) 0-homologous oriented 
subgraph, at least one edge of Ti corresponds to a stem of U. Let s be the last stem of 
U corresponding to an edge of Ti that is reattached by the complete closure procedure. 
Consider the step where s is reattached. As the root angle (and thus /o) is in the special 
face (see the terminology of Section]^, the special face is in the region defined by Xq. 
Thus it is on the left of s when it is reattached. This contradicts the fact that U is 
balanced. Thus D is the minimal HTC Schnyder wood w.r.t. /q. □ 


10 The lattice of HTC Schnyder woods 

In this section, we push further the study of HTC Schnyder woods in order to remove 
the root and the balanced property of the unicellular maps considered in Theorem and 
obtain a simplified bijection in Theorem of Section [m 
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Consider a toroidal triangulation G given with a crossing Schnyder wood. Let Dq 
be the corresponding 3-orientation of G. Let /o be any face of G. Recall from Section 
that 0{G) denotes the set of all the orientations of G that are homologous to Dq. The 
elements of 0{G) are the HTC Schnyder woods of G and {0(G), <fQ) is a distributive 
lattice. 

We need to reduce the graph G. We call an edge of G rigid w.r.t. 0(G) if it has 
the same orientation in all the elements of 0(G). Rigid edges do not play a role for the 
structure of 0(G). We delete them from G and call the obtained embedded graph G. 
Note that this graph is embedded but it is not necessarily a map, as some faces may not 
be homeomorphic to open disks. Note also that G might be empty if all the edges are 
rigid, i.e. 10(^)1 = 1 and G has no edge but a unique face that is all the surface. 


Lemma 11 Given an edge e of G, the following are equivalent: 

1. e is non-rigid 

2. e is contained in a ^-homologous oriented subgraph of Do 

3. e is contained in a 0-homologous oriented subgraph of any element of 0(G) 


Proof. (1 3) Let D G 0(G). If e is non-rigid, then it has a different orientation in 

two elements D', D" of 0(G). Then we can assume by symmetry that e has a different 
orientation in D and D' (otherwise in D and D” by symmetry). Since D, D' are homol¬ 
ogous to Do, they are also homologous to each other. So T = Z) \ is a 0-homologous 
oriented subgraph of D that contains e. 

(3 =► 2) Trivial since Do G 0(G) 

(2 1) If an edge e is contained in a 0-homologous oriented subgraph T of Dq. 

Then let D be the element of 0(G) such that T = Do\D. Clearly e is oriented differently 
in D and Dq, thus it is non-rigid. □ 


By Lemma 11, one can build G by keeping only the edges that are contained in a 


0-homologous oriented subgraph of Do. Note that this implies that all the edges of G 
are incident to two distinct faces of G. Denote by F the set of oriented subgraphs of 
G corresponding to the boundaries of faces of G considered counterclockwise. Let fo be 
the face of G containing fo and Fo be the element of F corresponding to the boundary 
of fo- Let F' = F\Fo. The proof of [121 Theorem 7] shows that the elements of F' are 
sufficient to generate the entire lattice (0(G), </o), i-e. two elements D, D' of 0(G) are 
linked in the Hasse diagram of the lattice, with D <f^ D', if and only if D \ D' G F'. 


Lemma 12 For every element F ^ F there exists D in 0(G) such that F is an oriented 
subgraph of D. 

Proof. Let F € F. Let D be an element of 0(G) that maximize the number of edges 
of F that have the same orientation in F and D (i.e. that maximize the number of 
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edges of D oriented counterclockwise on the border of the face of G corresponding to 
F). Suppose by contradiction that there is an edge e of F that does not have the same 
orientation in F and D. Edge e is in G so it is non-rigid. Let D' G 0{G) such that 
e is oriented differently in D and D'. Let T = D \D'. By [T21 Claim 1 of the proof 
of Theorem 7], there exists edge-disjoint oriented subgraphs Ti,... ,Tk of D such that 
4>{T) = for 1 < i < /c, there exists Xi C F' and e* G { — 1,1} such 

that 4>{Ti) = ei W.l.o.g., we can assume that e is an edge of Ti. Let D” 

be the element of 0{G) such that Ti = D \ D”. The oriented subgraph Ti intersects F 
only on edges of D oriented clockwise on the border of F. So D" contains strictly more 
edges oriented counterclockwise on the border of the face F than D, a contradiction. 
So all the edges of F have the same orientation in D. So F is a 0-homologous oriented 
subgraph of F. □ 


By Lemma 12, for every element F € F' there exists D in 0{G) such that F is an 


oriented subgraph of D. Thus there exists D' such that F = D\D' and D, D' are linked 
in the Basse diagram of the lattice. Thus the elements of F' form a minimal set that 
generates the lattice. 

Let Fmax (resp. Fmin) be the maximal (resp. minimal) element of {0(G), </o)- 


Lemma 13 Fq (resp. —Fq) is an oriented subgraph of (resp. D^\a). 


Proof. By Lemma 
D. Let T = F \ F 
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there exists F in 0{G) such that F is an oriented subgraph of 
. Since D <fg Fmax, the proof of [l2l Theorem 7], shows that 
the characteristic flow of T can be written as a combination with positive coefficients of 
characteristic fiows of F', i.e. 0(r) = YIfcF ^f4>{F) with A G L So T is disjoint 


from Fq. Thus Fq is an oriented subgraph of F^ax- The proof is similar for F^ 


□ 


Note that the above three lemmas hold in a more general context than just 0(G). 
Actually they hold for any lattice of homologous orientations on an oriented surface 
(see [H]). From now on we use some specific properties of the object considered in this 
paper, i.e. HTC Schnyder woods. 

Lemma 14 Consider an orientation D in 0(G) and a closed walk W of G. If on the 
left (resp. right) side of W, there is no incident edges of G, and no outgoing incident 
edges of D, then W is a contractible triangle with its contractible region on its left (resp. 
right) side. 

Proof. Consider a closed walk W of G such that on its left side there is no incident edge 
of G, and no outgoing incident edges of F. Let k be the length of W. Let Wi^ft be the 
edges of F that are incident to the left side of W. By assumption they are all entering 
W. Note that W cannot cross itself otherwise it has at least one incident edge of G on 
its left side. However it may have repeated vertices but in that case it intersects itself 
tangencially on the right side. 
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Suppose first that W is non-contractible. Then consider the closed walk W* of the 
dual orientation D* that is obtained by considering all the dual edges of VF/e/t with their 
corresponding orientation. Since all the edges of Wieft are entering W we have that W* 
is an oriented closed walk. Moreover it is non-contractible and thus contains an oriented 
non-contractible cycle, a contradiction to Lemma So W is contractible. Since W can 
intersect itself only tangencially on the right side, the region delimited by W and located 
on its left side is connected. 

Suppose that W has its contractible region on its left side. Consider the graph G' 
obtained from G by keeping only the vertices and edges that lie in the contractible region 
delimited by W, including W. The vertices of W appearing several times are duplicated 
so that G' is a plane triangulation of a A:-cycle. Let n',m', f be the number of vertices, 
edges and faces of G'^ By Euler’s formula, n' — m! + f = 2. All the inner faces have 
size 3 and the outer face has size k, so 2m' = 3(/' — 1) -|- A:. All the inner vertices 
have outdegree 3 as we are considering a 3-orientation of G. All the edges of Wi^ft are 
oriented toward W, and there are k outer edges, so m' = 2>{n' — k) + k. Combining these 
three equalities gives k = 3, i.e. VE is a triangle and the lemma holds. 

Suppose now that W has its contractible region on its right side. Then similarly as 
above, consider the graph G' obtained from G by keeping all the vertices and edges that 
lie in the contractible region delimited by W, including W. This time the vertices of W 
appearing several times are not duplicated. Since W can intersect itself only tangencially 
on the right side, we have that G' is a plane map whose outer face boundary is W and 
whose interior is triangulated. As above, let n',m',f' be the number of vertices, edges 
and faces of G'. By Euler’s formula, n' — m' + f = 2. All the inner faces have size 3 
and the outer face has size k, so 2m' = 3{f' — 1) -|- A:. Since there is no outgoing incident 
edges of D on the left side of FT, all the vertices of G' have outdegree 3 and m' = 3n'. 
Combining these three equalities gives k = —3, a contradiction. □ 


The boundary of a face of G may be composed of several closed walks. Let us 
call quasi-contractible the faces of G that are homeomorphic to a disk or to a disk 
with punctures. Note that such a face may have several boundaries (if there is some 
punctures) and then the face is not contractible, but exactly one of these boundaries 
contains all the other in its contractible region. Let us call outer facial walk this special 
boundary. Then we have the following: 


Lemma 15 All the faces of G are quasi-contractible and their outer facial walk is a 
(contractible) triangle. 


Proof. Suppose by contradiction that there is a face / of G that is not quasi-contractible 
or whose outer facial walk is not a contractible triangle. Let F be the element of F 

there exists an orientation D in 
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corresponding to the boundary of /. By Lemma 
0{G) such that F is an oriented subgraph of D. 

All the faces of G, are contractible triangles. Thus / is not a face of G and contains in 
its interior at least one edge of G. Start from any such edge e and consider the left-walk 
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W = (ei)i>o of D obtained by the following: if the edge e* is entering a vertex v, then 
Cj+i is choosen among the three edges leaving v as the edge that is on the left coming 
from ej (i.e. the first one while going clockwise around v). Suppose that for f > 0, edge 
Cj is entering a vertex v that is on the border of /. Recall that by definition F is oriented 
counterclockwise according to its interior, so either Cj+i is in the interior of / or e^+i is 
on the border of /. Thus W cannot leave /. 

Since G has a finite number of edges, some edges are used several times in W. 
Consider a minimal subsequence W' = ,..., such that no edge appears twice and 

Cfc = e^+i- Thus W ends periodically on the sequence of edges e^,..., e^. By Lemma 
all the closed walks that are part of F have some outgoing incident edges of D on their 
left side. Thus we have that W contains at least one edge that is not an edge of F, thus 
it contains at least one rigid edge. 

By construction, all the edges on the left side of W are entering. Suppose that W 
is not contractible. Then the oriented closed walk of the dual orientation D* that is 
obtained by considering all the dual edges of its incident edges on the left side gives 
a contradiction to Lemma So W' is contractible. So it is a 0-homologous oriented 
subgraph of D, thus all its edges are non-rigid by Lemma HH a contradiction. □ 
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By Lemma 15 every face of G is quasi-contractible and its outer facial walk is a 


contractible triangle. So G contains all the contractible triangles of G whose interiors 
are maximal by inclusion, i.e. it contains all the edges that are not in the interior of a 
separating triangle. In particular, G is non-empty and |0(G')| > 2. The status (rigid or 
not) of an edge lying inside a separating triangle is determined as in the planar case: such 
an edge is rigid if and only if it is in the interior of a separating triangle and incident to 
this triangle. Thus an edge of G is rigid if and only if it is in the interior of a separating 
triangle and incident to this triangle. 

Since (0 (G),</q) is a distributive lattice, any element D of 0{G) that is distinct 
from Umax and Umin contains at least one neighbor above and at least one neighbor 
below in the Basse diagram of the lattice. Thus it has at least one face of G oriented 
counterclockwise and at least one face of G oriented clockwise. Thus by Lemma 
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it contains at least one contractible triangle oriented counterclockwise and at least one 
contractible triangle oriented clockwise. Next lemma shows that this property is also 
true for Umax and T>min- 


Lemma 16 In Umax (resp. UminJ there is a counterclockwise (resp. clockwise) con¬ 
tractible triangle containing /o, and a clockwise (resp. counterclockwise) contractible 
triangle not containing /q. 


Proof. By Lemma |15[ fn is quasi-contractible and its outer facial walk is a contractible 
triangle T. By lemma [l^ Fq is an oriented subgraph of Umax- Thus T is oriented 
counterclockwise and contains /q. The second part of the lemma is clear since |0(G)| > 2 
so Umax has at least one neighbor below in the Basse diagram of the lattice. Similarly 
for Dmin- n 
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Thus by above remarks and Lemma 16, all the HTC Schnyder woods have at least 
one triangle oriented counterclockwise and at least one triangle oriented clockwise. Note 
that this property does not characterize HTC Schnyder woods. Figurej^gives an example 
of a Schnyder wood that is not HTC but satisfies the property. Note also that not all 
Schnyder woods satisfy the property. The right of Figure is an example of a Schnyder 
wood that is no HTC and has no oriented contractible triangle. 

Lemma [T6| is used in the next section to obtained a bijection with unrooted unicellular 
maps. 


11 Bijection with unrooted unicellular maps 


To remove the root and the balanced property of the unicellular maps considered in 
Theorem we have to root the toroidal triangulation more precisely than before. We 
say that an angle is not in the clockwise interior of a separating triangle if it is not in 
its contractible region, or if it is incident to a vertex v of the triangle and situated just 
before an edge of the triangle in counterclockwise order around v (see Figure 19). 



Figure 19: Angles that are in a separating triangle but not in its clockwise interior. 

Consider a toroidal triangulation G. Consider a root angle uq that is not in the 
clockwise interior of a separating triangle. Note that the choice of oq is equivalent to the 
choice of a root vertex vq and a root edge cq incident to vq such that none is in the interior 
of a separating triangle. Consider the orientation of the edges of G corresponding to the 
minimal HTC Schnyder wood w.r.t. the root face /q. By Lemma [161 there is a clockwise 
triangle containing /q. Thus by the choice of oq, the edge cq is leaving the root vertex vq. 
This is the essential property used in this section. Consider the output U of Algorithm 
PS on (G, oo). Since cq is leaving vq and oq is just before cq in counterclockwise order 
around vq, the execution of Algorithm PS starts by Case 2 and eo corresponds in U 
to a stem sq attached to vq. We call this stem sq the root stem. 

The recovering method defined in Theorem says that sq is the last stem reattached 
by the procedure. So there exists a sequence of admissible triples of U (see the termi¬ 
nology and notations of Section such that sq belongs to the last admissible triple. 
Let Uq = U and for 1 < A: < 2n — 2, the map is obtained from Uk-i by closing any 
admissible triple that does no contain sq. As noted in Section[^ the special face of U 2 n -2 
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is a quadrangle with exactly one stem. This stem being sq, we are in the situation of 
Figure!^ 



Figure 20: The situation just before the last stem (i.e. the root stem) is reattached 

Consequently, if one removes the root stem sq from U to obtain an unicellular map 
U' with n vertices, n + 1 edges and 2n — 2 stems, one can recover the graph U 2 n -2 by 
applying a complete closure procedure on U' (see example of Figure [^. Note that then, 
there are four different ways to finish the closure of U 2 n -2 to obtain an oriented toroidal 
triangulation. This four cases corresponds to the four ways to place the (removed) root 
stem in a quadrangle, they are obtained by pivoting Figureby 0°, 90°, 180° and 270°. 
Note that only one of this four cases leads to the original rooted triangulation G, except 
if there are some symmetries (like in the example of Figure [2T|). 


2 6 3 2 6 3 



Non-rooted output Partial recovering 


Figure 21: Example of Kj where the root angle, the root stem and the orientation 
w.r.t. the root angle have been removed from the output of Figure The complete 
closure procedure leads to a quadrangular face. 

Let U{n) denote the set of (non-rooted) toroidal unicellular maps, with exactly n 
vertices, n -|- 1 edges and 2n — 2 stems satisfying the following: every vertex has exactly 
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2 stems, except if the map is hexagonal, the two corners having exactly 1 stem each, 
and if the map is square, the only corner having no stem at all. Note that the output 
of Theorem]^ on an appropriately rooted toroidal triangulation is an element of hl{n) 
when the root stem is removed. 


Note that an element U' of lA{n) is non-rooted so we cannot orient automatically its 
edges w.r.t. the root angle like in Section]^ Nevertheless one can still orient all the 
stems as outgoing and compute 7 on the cycles of U' by considering only its stems in 
the counting (and not the edges nor the root stem anymore). We say that an unicellular 
map of U{n) satisfies the 70 property if 7 equals zero on its (non-contractible) cycles. 
Let us call (n) the set of elements of U (n) satisfying the 70 property. 

A surprising property is that an element U' of hl{n) satisfies the 70 property if and 
only if any element U of lAr (n) obtained from U' by adding a root stem anywhere in U' 
satisfies the 70 property (note that in U we count the edges and the root stem to compute 
7 ). One can see this by considering the unicellular map of Figure 22 It represents the 


general case of the underlying rooted hexagon of U. The edges represent in fact paths 
(some of which can be of length zero). One can check that it satisfies 7 equals zero on 
its (non-contractible) cycles. It corresponds exactly to the set of edges that are taken 
into consideration when computing 7 on [/ but not when computing 7 on U'. Thus it 
does not affect the counting (the tree-like parts are not represented since they do not 
affect the value 7 ). So the output of Theorem on an appropriately rooted toroidal 
triangulation is an element of hly^{n) when the root stem is removed. 



Figure 22: The parts of the unicellular map showing the correspondence while computing 
7 with or without the orientation w.r.t. the root plus the root stem. 

For the particular case of Kt, the difference between the rooted output of Figure 
and the non-rooted output of Figure [M] is represented on Figure (one can superimpose 
the last two to obtain the first). One can check that these three unicellular maps (rooted, 
non-rooted and the difference) all satisfy 7 equals zero on their cycles. 

There is an “almost” four-to-one correspondence between toroidal triangulations on 
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Figure 23: The difference between the rooted output of Figure and the non-rooted 
output of Figure 

n vertices, given with a root angle that is not in the clockwise interior of a separating 
triangle, and elements The “almost” means that if the automorphism group 

of an element U of Ll.yg{n) is not trivial, some of the four ways to add a root stem 
in U are isomorphic and lead to the same rooted triangulation. In the example of 
Figure one can root in fonr ways the quadrangle but this gives only two different 
rooted triangulations (because of the symmetries of K^). We face this problem by 
defining another class for which we can formulate a bijection. 

Let T (n) be the set of toroidal maps on n vertices, where all the faces are triangles, 
except one that is a quadrangle and which is not in a separating triangle. Then we have 
the following bijection: 

Theorem 7 There is a bijection between T{n) and 

Proof. Let a (for “add”) be an arbitrarily chosen mapping defined on the maps G' of 
T(n) that adds a diagonal eo in the quadrangle of G' and roots the obtained toroidal 
triangnlation G at a vertex vq incident to cq (this defines the root angle ag situated 
just before eg in counterclockwise order around vg). Note that the added edge cannot 
create a separating 2 -cycle, since otherwise the quadrangle wonld be in a separating 
triangle. Moreover the root angle of G is not in the clockwise interior of a separating 
triangle. Thus the image of a is in Tf{n), the snbset of Tr{n) corresponding to toroidal 
triangnlations rooted at an angle that is not in the clockwise interior of a separating 
triangle. 

Let 2/' ^ (n) be the elements of ^r, 6 , 7 o (^) have their root angle jnst before a 

stem in counterclockwise order aronnd the root vertex. Consider the mapping g, defined 
in the proof of Theorem By above remarks and Theorem the image of g restricted 
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to Tr{n) is in ^r 6 7 o(^)' “remove”) be the mapping that associates to an 

element of ^ (n) an element of (n) obtained by removing the root angle and its 
corresponding stem. Finally, let h = r o g o a which associates to an element of 7~(n) an 
element of U^Q{n). Let us show that /i is a bijection. 

Consider an element G' of T (n) and its image U’ by h. The complete closure proce¬ 
dure on U' gives G' thus the mapping h is injective. 

Conversely, consider an element U' of U-y^ (n). Apply the complete closure procedure 
on U'. At the end of this procedure, the special face is a quadrangle whose angles are 
denoted ..., a^. We denote also by ..., the corresponding angles of U'. For 
i G {1,..., 4}, let be the element of Ur{n) obtained by adding a root stem and a root 
angle in the angle a* of U\ with the root angle just before the stem in counterclockwise 
order around the root vertex. Note that by the choice of a*, the C/* are all balanced. By 
above remarks they also satisfy the 70 property and thus they are in U'^ ^ (n). 

By the proof of Theorem the complete closure procedure on gives a triangula¬ 
tion G* of Tr{n) that is rooted from an angle Oq not in the strict interior of a separating 
triangle and oriented according to the minimal HTC Schnyder wood w.r.t. the root face. 
Moreover the output of Algorithm PS applied on (G®, Og) is W. Since in W, the root 
stem is present just after the root angle, the first edge seen by the execution of Algo¬ 
rithm PS on (G^jOg) is outgoing. So oq is not in the clockwise interior of a separating 
triangle (in a 3-orientation, all the edges that are in the interior of a separating triangle 
and incident to the triangle are entering the triangle). Thus the G® are appropriately 
rooted and are elements of T^{n). Removing the root edge of any Gj, gives the same 
map G' of T{n). Exactly one of the Gi is the image of G' by the mapping a. Thus the 
image of G' by h is U' and the mapping h is surjective. □ 

A nice aspect of Theorem [^comparing to Theorem is that the unicellular maps that 
are considered are much simpler. They have no root nor balanced property anymore. 
It would be great to use Theorem to count and sample toroidal triangulations. The 
main issue comparing to the planar case is the 70 property. 


12 Higher genus 

The key lemmas that make the encoding method presented in this paper work are Lem¬ 
mas!^ and Note that Theorem]^ is proved in a very general form in |12] . So one can 
consider a minimal element in the lattice and get the equivalent of Lemma Things are 
more complicated for Lemma since the existence of Schnyder woods in higher genus 
is only conjectured when g > 2 and moreover we have no idea of what would be the 
generalization of crossing property and thus HTC Schnyder woods. Nevertheless one 
can hope to find orientations satisfying the conclusion of Lemma and thus apply the 
same encoding method as here. This is what we discuss below. 

Recently, Albar, the second author and Knauer [1] proved the following: 
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Theorem 8 (|T]) A simple triangulation on a genus g > 1 orientable surface admits 
an orientation of its edges such that every vertex has outdegree at least 3, and divisible 
by 3. 

Theorem is proved for simple triangulation but we believe it to be true for all 
triangulations. Moreover we hope for a possible generalization satisfying the conclusion 
of Lemma [2l 

Conjecture 1 A triangulation on a genus g > 1 orientable surface admits an orienta¬ 
tion of its edges such that every vertex has outdegree at least 3, divisible by 3, and such 
that there is no oriented non-contractible cycle in the dual orientation. 

Even if Conjecture is true, more efforts should be made to obtain a bijection since 
there might be several minimal element satisfying the conjecture and one has to identify 
a particular one (like the minimal HTC Schnyder wood in our situation). 
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